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(N ■ Abstract 

We study a general car-following model of traffic flow on an infinitely long single-lane road, which 
assumes that a car's acceleration depends on its own speed, the headway between it and the car 
ahead, and the rate of change of headway, but makes minimal assumptions about the functional 

(N : 

£Nj , form of that dependence. The velocity of uniform steady flow, as a function of headway, can be 

found implicitly from the acceleration function, and the linear stability criterion for that flow can be 
expressed simply in terms it. Crucially, and unlike in previously analyzed car-following models, the 
threshold of absolute stability does not generally coincide with an inflection point in the steady- 
state velocity function. The Burgers and Korteweg-deVries equations can be derived under the 
usual assumptions, but the modified Korteweg-deVries equation arises only when the threshold of 
■ absolute stability does coincide with an inflection point. Otherwise, the Korteweg-deVries equation 



O 

(N 



\f} ', continues to apply near absolute stability, while near the inflection point one obtains the modified 

Korteweg-deVries equation plus an extra, quadratic term. Corrections to the Korteweg-deVries 
equation "select" a single member of the a one-parameter set of one-soliton solutions by driving 
a slow evolution of the soliton parameter either away from or toward that member. In previous 
models this selected soliton has always marked the threshold of a finite-amplitude instability of 
linearly stable steady flow, but here it can alternatively be a stable, small-amplitude jam that 
occurs when steady flow is linearly unstable. The latter is the case when the Korteweg-deVries 
equation applies near absolute stability. The new, augmented mKdV equation, which holds near an 
inflection point, admits a continuous family of kink solutions, like the modified Korteweg-deVries 
equation. We derive the selection criterion arising from the corrections to this equation. 

PACS numbers: 05.45.Yv, 45.70.Vn, 47.20.Ky, 47.54.-r 
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I. INTRODUCTION 



Car-following models constitute a class of highly idealized, "microscopic" models of ve- 
hicular traffic, which attempt to capture the characteristics of traffic flow by modeling the 
behavior of the individual cars. A typical car-following model describes a single line of cars 
traveling along a long, uniform, straight or circular road. Each driver is assumed to adjust 
his or her speed according to what the next car ahead is doing, and all drivers are taken to 
behave identically. It is well known that uniform, steady traffic flow can be unstable against 
the formation of slow- moving traffic jams; these models are used to identify conditions under 
which these jams form and to investigate their properties. 

An importantprototype of this class of models is the "optimal velocity model" (OVM) 
of Bando et al. [1], which is embodied in the equations 

r^ = V s (x n+1 -x n )-v n , v n = dx n /dt. 

Here the cars move in the positive x direction, they are numbered consecutively, with car 
n + 1 ahead of car n, and x n denotes the position of, say, the back of car n. The "optimal 
velocity" function V s (h) is the speed at which a driver prefers to drive when the next car is 
ahead by a distance, or "headway," of h. The model assumes that each driver relaxes his 
or her car's speed to that preferred velocity with some fixed time constant r. The optimal 
velocity function also acquires a second interpretation in this model, as it is trivial to show 
that V s {h) is the steady-state speed of uniform traffic. That is, it is the speed at which every 
car moves when the cars are spaced an equal distance h apart. 

A model like the OVM, of course, ignores many features of the behavior of real human 
drivers, who not only have different reaction times, but also have individual preferred velocity 
functions which can change because of road conditions, fatigue, distractions, and the like. 
However, such a model would give a better description of a line of cars controlled by on-board 
adaptive cruise control systems, which would use a small radar to monitor the headway to 
the next car. The equation of motion would then be implemented by the programming of 
the cruise control. 

The OVM was quickly generalized j^] to account for the fact that the proper response to 
make to a car that is a distance h ahead must also depend on that car's speed - the response 
of car n must be very different if the next car in line is stopped than it would be if that car 
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is moving at some speed near v n . A particularly simple model that incorporates this effect 
is the "full velocity difference model" (FVDM) of Jiang et al. j^], given by 

r ^ = V s (x n+1 - x n ) - v ri + A (v n+1 - v n ), 

where A is a constant which specifies the relative importance of the size of the headway 
and its rate of change in determining the response of car n. Note that in this model, the 
function V s (h) retains its interpretation as the steady-state speed of uniformly spaced traffic, 
but it is not the speed to which an individual driver tends in all circumstances. Instead, an 
individual driver tends to adjust his or her speed toward some weighted average of V s (h) and 
v n+ i, rather than just toward V s {h). Tian, Jia, and Li |4| suggested replacing the v n+ \ — v n 
in the FVDM with a general function of v n+ \ — v n , naming their model the "comprehensive 
optimal velocity model" (COVM); they carried out calculations in the specific case where 
the function was a hyperbolic tangent. 

A number of car-following models have been proposed in the literature; for reviews see 
5KZ[. Analyses of different car-following models yield a range of results that are strikingly 
similar from model to model: 

(1) Steady, uniform flow with a fixed headway h — A is possible for a continuous range 
of headway values. This flow is linearly stable or linearly unstable depending on whether 
V^A), the derivative of the preferred velocity function, is below or above some critical value 
which depends on the parameters of the model. (The preferred velocity function is often 
taken to have a sigmoidal, hyperbolic-tangent-like shape, in which case this means that 
light traffic, with a large spacing A, is stable, and very heavy traffic, with small A, is also 
stable albeit slow, while instability can happen for intermediate traffic densities. It is then 
possible to have absolute stability, when the slope V^'(A) never exceeds the critical value.) 
When V^'(A) first crosses the critical value, the first disturbances to steady flow that become 
unstable are those with vanishingly small wave numbers. 

(2) If one assumes that the headway between cars varies only slightly from a uniform 
value A, and that it varies slowly from car to car, then one can derive a Burgers equation 
for headway as a function of time and car number, as has been done for the OVM [8j, the 
FVDM [9[, and the COVM 4|. For those models, the coefficient of the nonlinear term 
in the Burgers equation is given by V"(A), and the diffusion coefficient is proportional to 
the difference between V^'(A) and its critical value. When uniform traffic is linearly stable, 



the diffusion coefficient is positive, and the Burgers equation then describes how slowly- 
varying deviations from uniformity diffuse away (via an interesting intermediate-asymptotic 



evolution [10|). When uniform flow is linearly unstable, the diffusion coefficient is negative, so 
the Burgers equation becomes a backwards diffusion equation. Smooth initial perturbations 
then grow sharper, rapidly invalidating the assumptions under which the Burgers equation 
is derived. 

(3) When V^'(A) is close to its critical value, one can reduce the model to a Korteweg- 
deVries (Kd V) eq uation for the headway, plus small corrections; again, this has been done 



for the OVM 



the FVDM (StllJ, and the COVM [4|. The KdV equation is completely 
integrable, and has soliton (and multi-soliton) solutions that can be written in closed form. 
The correction terms, which break the integrability of the equation, have a very suggestive 
effect on the one-soliton solutions: depending on the signs of certain coefficients in the 
correction terms, there may be at most one single member of the continuous family of one- 
soliton solutions that survives as a localized, shape-preserving solution. Other members of 
the one-soliton family are slowly driven either toward or away from this one member jl^ . \v\ , 
again depending on signs of correction coefficients. In general, however, there is a surviving 
soliton solution only when steady flow is linearly stable, and it marks the threshold of a 
finite-amplitude instability of that flow. 

(4) Near the onset of absolute stability - that is, when there is at most a narrow range of 



parameters for which linear instability can occur - the headway is described 



Dy a modified 



Korteweg-deVries (mKdV) equation, plus correction terms, in the OVM 15], the FVDM 
9), Q, and the COVM [4]. This equation is also completely integrable; instead of localized 
soliton solutions, however, it has kink (and multi-kink) solutions which look like the transi- 
tion zones between regions of different (linearly stable) traffic densities. As is the case for 
one-soliton solutions of the KdV equation, there is a continuous family of one-kink solutions 
of the mKdV equation, and the correction terms select a single member of this family. 

The purpose of this paper is to test the generality of these results. For instance, the 
reduced equations (Burgers, KdV, mKdV) are derived from fairly long perturbation series 
about the steady-flow state. Are the results of these derivations dependent on the fact 
that the OVM and FVDM assume that the acceleration of linear function of its 

velocity? In the FVDM, the parameter A, which specifies the relative importance of the size 
of the gap and its rate of change, is taken to be a constant. Clearly it is unrealistic to expect 



the speed of car n + 1 to be relevant if it is, say, several kilometers ahead of car n. Do any 
important results change if it, or the time constant r, is allowed to depend on headway, as 
suggested by Gasser et al. [16|? To answer questions such as these, we analyze a general 
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possible car-following model, embodied in the equation 

^ = A(x n+1 - x n , v n+1 - v n ,v n ). (1.1) 

Here the acceleration function A(h, h, v) is a general function of the speed v = v n of the 
car under consideration, the headway h = x n+ \ — x n between it and the next car ahead, 
and the rate h = v n+ \ — v n at which the headway is changing. To be a realistic model of 
driver behavior, this function must satisfy some minimal conditions which we will discuss 
in Section [III Note that the model still assumes a single line of cars, so that the behavior 
of car n depends only on the single car ahead of it, and it also assumes a uniform road, 
in that there is no position dependence in the acceleration function. The model does not 
include explicit, finite time delays or look-ahead or look-back effects (viz., a dependence of 
the behavior of car n on car n + 2 or car n — 1 ) . 

We find that many, but not all, of the results found in specific car-following models to date 
are in fact universal at the level of generality of ( 11. ip . Specifically, the linear stability results 
and the derivations of the Burgers and KdV equations come out as usual. However, a new 
term appears in the corrections to the KdV equation, and this term opens up the possibility 
of new behaviors not seen in its absence. The same term shows up in the leading order in 
the derivation that usually produces the mKdV equation, so that the mKdV equation is 
not generic. In the mKdV context, the possibility of having an extra term of this form was 



anticipated by Komatsu and Sasa 151 ]. and it was also found in a model including look-back 



id- 



by Hayakawa and Nakanishi 

In Section [TT] we discuss the conditions that the function A(h, h,v) must satisfy in order 
for (11.11) to be a reasonable model of traffic flow, and we carry out the linear stability analysis 
of steady traffic flow with a uniform headway A. We show how the analysis can be phrased 
in terms of the steady-state velocity function V S (A), which, although it does not appear 
in the formulation of the model, can be defined implicitly from the acceleration function 
A(h,h,v). We find that the linear stability results quoted above continue to hold at the 
level of generality of ( II. ip . but that the condition for absolute stability must be stated in 
a more general way than usual. In Section [TTT] we expand (II .ip about the uniform steady 



state, assuming that the headway varies only slowly from car to car. From the leading order 
of this expansion we find that the Burgers equation appears exactly as described above. In 
Section HVl we modify the expansion from Section HTT1 to apply to a uniform steady state that 
is close to the onset of instability. We find that the KdV equation appears at leading order 
as described above, but that an extra term appears in the coefficient of one of the corrections 
to the KdV equation. We then carry out a multiple-time-scales analysis of the one-soliton 
solutions of the leading-order KdV equation. We find that when the new term is absent, the 
unique one-soliton solution that is preserved by the corrections marks the onset of a finite- 
amplitude instability of a linearly stable uniform flow; with the new term present, however, 
it is possible instead for this one-soliton solution to represent a small- amplitude traffic jam 
that is stable when uniform flow is slightly unstable. In Section |V] we again modify the 
expansion from Section IHH this time to apply to a uniform steady state that is close to 
an inflection point of the steady-state velocity function. This normally corresponds to the 
onset of absolute stability, but at the level of generality of ( 11. II) this need not be the case. 
We find now that the new term appears at leading order - that is, the flow is described by 
an evolution equation that is not the mKdV equation. (The new term, however, does vanish 
if the inflection point coincides with the threshold of absolute stability.) Like the mKdV 
equation, this new equation also has a one-parameter family of kink solutions. Moreover, 
we show that it is still possible to carry out the selection analysis for these kink solutions, 
even though we have not been able to show that the new equation is integrable. We discuss 
our results in Section [VI] 



II. STEADY STATES AND LINEAR STABILITY 



Before deriving and analyzing the steady states of the model embodied in f ll.ip . we first 
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set forth some conditions which a realistic acceleration function A(h } h,v) must satisfy 
Here and throughout this paper we will assume that A(h,h,v) is differentiable as many 
times as necessary. For a given headway h and a given rate of change of headway h, a 
driver will be more prone to decelerate or less prone to accelerate the faster his or her car is 
traveling. Similarly, other things being equal a driver is more prone to accelerate the larger 



the headway is, or the more rapidly it is increasing. Thus we expect to have 

dA dA dA , s 

Tp^O, — >0, tt>0. 2.1 
ov dh dh 

We will in fact assume slightly more, namely that dA/dv is strictly negative, or at least that 

there are no finite ranges of v over which A(h, 0, v) is constant. 

To find simple steady-state solutions of the model, we assume that the cars are equally 

spaced, and that they all cruise at the same speed, so they remain equally spaced. If the 

headway in front of each car is A, then this uniform flow is a steady state provided the 

common speed of the cars is V S (A), which is defined implicitly by 

A(A,0,V s (A)) = 0. (2.2) 

From the first inequality in (12. ip we see that A(A,0,v) is a nonincreasing function of v; 

given our assumption above that A(h, 0, v) does not remain constant over any finite range of 

v, the steady flow speed for a given A will then be unique. Thus we have a one-parameter 

family of steady states of the model, labeled by the spacing A; of course, these steady-flow 

solutions may be stable or unstable. 

By differentiating (12.21) repeatedly with respect to A, we can obtain expressions for the 

derivatives of V s . For example, we find 

_ A(A,0,K(A)) 
VA) ~ A„(A,0,V s (A))' {2J} 

where, here and henceforth, primes denote derivatives with respect to A and subscripts 

denote partial derivatives, e.g. A^ = dA/dh. Given the restrictions (12. ip on the signs of the 

derivatives of A, we see that V^'(A) cannot be negative. By differentiating (I2.3p . we find 

v»tM - d " Vs - A hh + 2A hv v:(A) + A w V a a (A) 

Vs lAj " rfA 2 " A v ' 1 j 

in which (again, here and henceforth) is it to be understood that the partial derivatives are 

evaluated in the steady state {h — A, h — 0, v — V S (A)). In like manner we can also derive 

expressions for higher derivatives of V s . 

The linear stability of these steady-flow solutions has been analyzed at this level of gen- 



erality by Wilson 17| and by Orosz et al. [19|; we recapitulate that calculation here, mainly 



to fix notation and to make contact with the linear stability results for specific models in 
the literature. We proceed in the usual way, first writing 

x n = nA + V s {A)t + x n {t), (2.5) 



where x n represents an infinitesimal deviation of the car positions from the steady-flow 
solution. We substitute this into the equation of motion (ll.ip and linearize in the x n to get 

d 2 x n _ V^(A) I dx n \ d 

,,9 — \ x n+l ~ x n) 77 I rrl^n+l — x n)i l^-"J 

at z t t at t at 

where the new parameters r and A are defined by 

- = -A V (A, 0, V 3 (A)), - = A k (A, 0, V S (A)). (2.7) 
r r 



These definitions agree with the notation of Ou et al. (9], but unlike in that paper A and 
t can now be nontrivial functions of the uniform spacing A. We see from the inequalities 
(12. ip that they must be positive (or at least non- negative). Note for future reference that 
(12. 4p . with this notation, reduces to 

K"(A) = r(A)[A hh + 2A h X(&) + A VV V®{A)]. (2.8) 

Similarly, we find 

(7) = -14* + A»K(A)], 0)' = + ^K'(A)], (2.9) 

and so forth. 

Since the coefficients in the linearized evolution equation ( 12. 6 p are independent of car 
index n and time t, we write the deviation x n (t) as a linear combination of Fourier modes 
exp(ikn), with amplitudes that grow or decay exponentially in time, as exp[a(k)t\. We then 
find that the (complex) linear growth rates cr(k) are given by 

a 2 r(A) = V^(A)(e lk - 1) - a + \(A)a(e ik - 1). (2.10) 

A mode of wave number k is marginally unstable when its growth rate cr(k) is pure imaginary. 
Separating the real and imaginary parts of (I2.10p then shows, after some algebra, that this 
occurs for 

V:=^ [l + (l + 2A)tan 2 |], (2.11) 

and that when the parameters satisfy this condition, the two solutions of the quadratic 
(|2~T0|) are 

or = i(l + 2A)tan^ and or = - ( 1 + 2A sin 2 - J fl+itan-). (2.12) 
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Clearly, the second solution has negative real part, and so represents a stable mode; the first 
can easily be shown to acquire a positive real part when V s ' is larger than the critical value 
given in (12. lip . The steady state is unstable if any perturbation grows, so it is unstable if 
(for a fixed A) V s ' exceeds the smallest possible value of the right side of (12. lip , which is 
reached at k = 0. Thus we find that the steady state is 

linearly unstable for V S \A) > fi c (A) = 1 + , (2.13) 

2r(A) 

and if V s ' is increased past the instability limit Q c , then it is the long-wavelength (k = 0) 
perturbations that first become unstable. 

Depending on the details of the acceleration function, it is possible for the steady-flow 
states to be "absolutely stable," in the sense that there is no A for which the flow is linearly 
unstable. This is the case if V S '(A) is always less than Q C (A). The threshold of absolute 
stability would be where only a single A is marginally stable, i.e. where V s — Q c has a 
maximum at a value of zero. For this to happen, there must be a A for which the value 
and slope of V s ' match the value and slope of Q c , and for which V"' is less than (or more 
negative than) the value of Q". In the OVM, the FVDM, and the COVM the parameters r 
and A are constants independent of A, and so Q c is a constant also. Thus in those models 
the threshold of absolute stability occurs at an inflection point of V s - but in more general 
models this need not be the case. 

Whether or not V s ' is near the threshold of instability, it is easy to extract the behavior of 
long- wavelength perturbations from (I2.10p . We immediately see that for small k, one of the 
solutions has a ~ — 1/r, so that this represents a decaying perturbation. The other solution 
can be found by expanding in powers of k, and is given by 



= V^k + 0(k 3 ), 

Wct = tV'-X-- 
2 



Vy + 0(k*). (2.14) 



From the real part, we see that this perturbation grows or decays according to whether V s ' is 
above or below the threshold given in (I2.13p . (When it is at the threshold, the neglected /c 4 
term turns out to be negative, so the perturbation then decays.) From the imaginary part, 
we see that the perturbation moves as it grows or decays, at a rate (in cars per unit time) of 
V^'(A) in the upstream direction relative to the cars. Since the cars are a distance A apart, 
at least before the perturbation develops, this is a distance per unit time of AV^'(A) relative 



to the cars. And since the cars are moving in the downstream direction at V S (A) relative 
to the road, the perturbation moves relative to the road at a velocity of approximately 
V S (A) — AVg(A), which may be positive (downstream) or negative (upstream). 

The stability criterion (I2.13P is identical to that obtained by Ou et al. |9[ for the FVDM, 
except for the more general definitions of V' a , r, and A given in (12.31) and (12.7ft . This should 
come as no surprise, since their model is essentially the general model expanded to linear 
order in the car velocities. Thus linear stability analysis is not sensitive to the generalizations 
implied in our evolution equation (11.11) . 

III. LONG- WAVE PERTURBATIONS 

As we saw from the linear stability analysis, when traffic conditions are such that steady 
flow is just on the verge of instability, it is the long-wavelength perturbations that are the 
first to grow. It is therefore of interest to approximate the full model in a way that focuses on 
these perturbations. To this end, we assume that the deviation of the cars' actual positions 
from their steady- flow positions varies slowly from one car to the next. From the small- A; 
linear stability result (I2.14p we see that each mode with a small wave number k grows or 
decays at a rate proportional to k 2 , while moving relative to traffic at a rate - in units of 
cars per unit time - of — V^'(A). Motivated by these considerations we write 

x n = nA + V s (A)t + f(z,T), (3.1) 

where 

z = e[n + Vt(A)t], T = e 2 t. (3.2) 

Here, e is a small but otherwise arbitrary parameter and / and its derivatives are taken 
to be of order unity. Thus in order for / to change by an amount of order unity, the car 
number n must change by an amount of order e _1 , so that e measures the slowness of the 
variation of the position deviation / along the line of cars. This ansatz amounts to assuming 
that there is a pattern of deviations of headway from the uniform value A, that this pattern 
propagates upstream through the line of cars at a rate V^'(A), and that its shape changes 
slowly as it propagates. 

If / is of order unity, then the velocities of the cars and the spacings between them deviate 
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from their steady-flow values by order e, 



v n = K(A) + eV;{A)f z + e 2 f T , (3.3) 
x n+1 - x n = A + ef z + ^e 2 f zz H , (3.4) 



where subscripts again denote partial derivatives. As we see, the function 

g(z,T) = f z (z,T) (3.5) 

then gives the leading-order deviation in both headway and velocity from the exact uniform- 
flow steady state; we will have occasion to express many of our results in terms of it. 

To derive an evolution equation for /, we substitute the ansatz (I3.1l) - (l3.2p into the equa- 
tion of motion (11.11) and expand in powers of e. The calculation is lengthy but straightfor- 
ward; carrying it out to fourth order yields (after canceling an overall factor e 2 ) 



h = Df zz + ~K'7f + e[C 11T f zT + C u f zzz + C 12T f z f T + C 12 f z f zz + C 13 f*] 

+ ^\p1\fzzzz + C 2 2afzfzzz + C 2 2bf zz + C 2 3f z fzz + C 24 f z H ] (3.6) 



The omitted terms in the e 2 correction all involve T-derivatives of /, and turn out not to 
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affect any of our analysis. The various coefficients are given by 

D=(Q C - VfrVt, 
Cut = A - 2tV s ', 



1 + 3A 

6 s ' 



C 12T = r(A hv + A V X) = -, 

T 



C\2 
Cl3 
C*21 
C*22a 



^[^w> + ( A hv + 2A hh )V: + 2A h X 2 ] = IK + 2tV:Q' c ), 

T -{A hhh + 3A hh X + 3A hv X 2 + A VVV V*) = ~ IK'" - —V" 

1 + 4A 

24 s ' 



\[A hh + (A hv + 3A hh )V: + 3A h X 2 ] = I 



V" + tV! 



1 + 3A 



T 



C 22b = ~(A hh + 4A h X + 4A h X 



(3.7a) 
(3.7b) 
(3.7c) 

(3.7d) 

(3.7e) 

(3.7f) 

(3-7g) 
(3.7h) 

(3.7i) 

(3-7j) 

v , (3.7k) 

where all derivatives of A are evaluated in the steady-flow state, and derivatives of V s , r, 
A, and Q c = (1 + 2A)/2r are evaluated at A. In the OVM and FVDM, the acceleration 
function A(h, h, v) is taken to be linear in the velocities h and v, with coefficients that are 
independent of h. For these models, then, r, A, and Q c are constants, and many of the 
coefficients above simplify considerably. As we will soon see, the fact that Cxi then reduces 
to V"/2 will prove particularly consequential. On the other hand, features that arise due to 
nonzero A-derivatives of r, A, and Q c will be novel. 

According to (13.61) . the position deviation / obeys a Burgers equation with corrections, 

1 



C 



23 



g[K" + 4rVX + tV:\A vv - AA hv + 4A hh )}, 

T ~[A hhh + 2(A hhv + A hhh X + (A hm + ±A hhv X 2 + 2A hm V: 3 } 



V 



3r' 



V" + 2XK + XV"(A VV - 2 A 



hv) 



C 



21 



T 

-^{Ahhhh 



±A hhh X + QA hhv X + AA hvv X + A vvvyVg ) 



1 

24 



At' 



V"" - —V'" + 6r - V" + 3tA vv V< 



f T = Df zz + -V:'f z +0(e), 



or, after differentiating with respect to z and using the definition (I3.5p . 

9t = Dg zz + V^gg z + 0(e), 



(3.8) 



(3.9) 
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This result has been obtained previously for the OVM the FVDM [9 : |, and the COVM 
J] ; the present calculation shows that the derivation of the Burgers equation, as well as the 
identification of the coefficients in terms of the steady velocity function V s , continue to be 
valid even at the level of generality of (II. ip . 

Near the onset of instability of steady traffic, D is small. Thus, in this regime, the 
neglected order-e corrections might have effects comparable in size to those of the linear 
term in the Burgers equation. We turn next to an examination of this situation. 



IV. NEAR THE ONSET OF INSTABILITY 



Suppose conditions are such that steady traffic flow is either slightly unstable or slightly 
stable, so that the parameter D in (I3.7ap is small. Let us then write 

V;(A) = n c (A) + e 2 6 } (4.1) 

where e is small and 5 is arbitrary; we could (but will not) choose 5 to be +1 when steady 
flow is slightly unstable or —1 when it is slightly stable, which would then give a specific 
meaning to the previously arbitrary small parameter e. From the small-fc linear stability 
result (I2.14p . we see that the unstable perturbations (for positive S) have wave numbers 
of order e. Then, once the ikV s in the imaginary part of the linear growth rate has been 
absorbed into the choice of reference frame in ( 13. 21) . the linear growth rate is then of order 
e 3 . Thus we retain the definition of z in (I3.2p . but replace the definition of T by 

T = eh. (4.2) 

Since V' a is near the onset of instability, we find from (I3.7al) that D is also of order e 2 . 
Thus in the Burgers equation (13.61) . either / must grow without bound or the fl term must 
eventually be balanced by the correction terms. In the latter case, the size of / must be 
only of order e, so that the quadratic term in the Burgers equation can be balanced by the 
linear term in the correction. (Other correction terms, being higher powers of /, would be 
of higher order.) Thus we modify the original ansatz ( 13. ip to read 

x n = nA + V s (A)t + ef(z,T), (4.3) 

We may now substitute this new ansatz into the evolution equation (11.11) and expand to 
fifth order in e. Equivalently, however, we can simply note that the new scalings have the 
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effect of multiplying each / in (13. 6p by an additional factor of e and also multiplying each 
T-derivative by an additional e. After canceling an overall factor of e 2 , this then gives, to 
order e, 



fT = C 11 f zzz + l -V'Jf 2 z +e 



1 + 2A 

-5fzz + C\yrfzT + Cyifzfzz + C21/22 



(4.4) 



2 

The leading order of this evolution equation is related to the Korteweg-deVries equation; if 
we take its z-derivative and rewrite it in terms of g = f z , then the leading terms take the 
standard KdV form 

1 + 2A, 



9t = C u g zzz + V"gg z + e 



— — Sg z + C UT g T + C 12 gg z + C 21 g zzz 



(4.5) 



-I z 



The correction terms can be simplified somewhat. Since V' s is close to Q c , we easily find 
Cut = _ (1 + A) to leading order. Replacing the gx in the correction with the leading-order 
terms, substituting the explicit expressions for the coefficients, and simplifying gives 

1 -I- 3A 1 -I- IX 

9t = —7—V' s g zzz + V:'gg z - ——e[A5g + (1 + 2A)lft„ + 2(1? - n' c )g% z . (4.6) 
o 8 

One-soliton solutions of the KdV equation - the leading order of (14. 6p - describe local 
increases or decreases in the car density that preserve their form while propagating along 
the line of cars. There is a one-parameter family of these solutions, given explicitly by 

g {0) (z,T; k) = Mk 2 sech 2 (kz + k\T), (4.7) 

with k being a free parameter, and 

M = 2[1 + 3A(A)]^^, u = \[l + 3A(A)]K(A). (4.8) 

Recall that z specifies the location of the soliton in a reference frame moving upstream at a 
rate (cars per unit time) of V^, so the effect of the u in this expression is to increase that rate 
by a factor 1 + 2(1 + 3A)e 2 fc 2 /3. The amplitude M is positive if V" is positive; in this case 
the soliton represents a local rarefaction of traffic, with both the headway and speed being 
higher in the center of the sech 2 lump. On the other hand, if V" is negative, then the soliton 
describes a local jam, a region where cars move more slowly and are closer together than 
average. For larger k the disturbance is larger in amplitude but smaller in spatial extent, and 
moves more rapidly through the line of traffic. The arbitrariness of k reflects the invariance 
of the KdV equation under the simultaneous rescalings z — > kz, T — > k 3 T, g — >■ k 2 g. This 
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invariance, in turn, reflects the arbitrariness of the small parameter e: replacing e with ke 
implements the rescalings. 

The correction term eg zz 5 in ( 14. 6 j) breaks the invariance of the equation under rescaling, 
and so one might expect that the corrections may break the continuous family of solutions 
(14 .7p down to a discrete set, or possibly a single solution. That is, we expect that only a 
discrete subset of the family of e = solutions actually represents the e — > limits of solutions 
of the full equation (14. 61) . Typically, one derives a solvability condition by linearizing ( 14.61) 



about the zero-order solution, and the surviving solutions have k values that satisfy that 
solvability condition. However, it is more instructive to calculate how the correction terms 
cause the value of k to evolve in time {13I . \h\ . The solvability condition gives no indication 
of whether the "selected" k is stable or unstable, while the multiple-time-scales approach 
does have something to say about this. 

To begin the calculation, we assume that we can write the first-order solution to (14. 6 p in 
the form 

g(z, T) = g*\z, T; k{eT)) + eg^(z, T). (4.9) 

That is, we allow the parameter k to evolve slowly in time. Substituting this into ( I4.6P and 
expanding to first order gives 

Cg {1) + = + (1 + 2\)V;g z °? + 2(K" - n' c )g^} zz , (4.10) 

where k denotes the derivative of k with respect to the slow time variable eT, and the linear 
operator £ is given by 

C = — - ±±*±V'— - V"-qW (4 11) 

dT 6 s dz* s dz 9 ■ [ ' 

The fact that there is a continuous family of solutions g( ) of the zero-order problem suggests 
that there will be a function that is annihilated by the adjoint operator C\ and in fact this 
function turns out to be g^ itself. Explicitly, with an inner product (-,-) that consists of 
integrating over all z and averaging over all T, the adjoint operator is 

t = _± l+3A y , q V)y»!L u 12) 

L &T + 6 V 'dz* +9 v 'dz- [ ] 

Applying this to g^ gives zero because g^ satisfies ( 14. 6 p to leading order. Thus if we take 
the inner product of g^ with both sides of (I4.10p . the term involving g^ vanishes and we 
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are left with 

(4.13) 

The usual solvability condition is simply the right hand side of this set equal to zero; it 
locates the k values which are not moved by the correction terms. Keeping the k term 
on the left reveals whether k's near those values are driven toward or away from them. 
Evaluating the integrals yields 

k = ^^^K {-511 + 2A(A), + 8,1 + 3A,A,]^p } 

(4.14) 

If we set Q' c = in this equation, then the k for which k vanishes agrees with the solvability 
result of Ou et al. |9| for the FVDM. 

Note that V" drops out of (I4.14p completely if Q' c vanishes. This is because for Q r c = 0, 
both the quadratic terms in (14. 6p . and none of the linear terms, carry factors of V". As a 
result, V" can be scaled out of the equation completely. For nonzero Q r c , on the other hand, 
there is a role for V" to play. 

In the OVM, the FVDM, and the COVM, fi c is a constant, so Vt' c = and the coefficient 
of k 5 on the right side of (14. 14j) is manifestly positive. Thus for 5 > - the range for which 
steady traffic flow is linearly unstable - (I4.14p drives the value of k higher and higher, in 
fact diverging in finite time. This means that the actual pattern of traffic flow evolves out 
of the regime in which the KdV equation (14. 6p is valid. On the other hand, for 5 < there 
is a nontrivial fixed point of (I4.14p with A; oc | (5 1 1//2 . \i k starts below this, then it decreases, 
decaying to zero as T _1//2 , while if k starts above the fixed point, then it is again driven to 
infinity. Thus the "selected" k actually marks the threshold of a finite- amplitude instability 
of steady traffic, even when steady flow is linearly stable. 

New possibilities arise when Q' c is nonzero, because the coefficient of k 5 in (14.141) can 
be negative. Should that be the case, there would be a nontrivial fixed point when 5 is 
positive. Thus when steady traffic flow is (slightly) linearly unstable, the instability leads 
to small-amplitude jams described by soliton solutions of the KdV equation (14. 6p . When 
steady flow is linearly stable, so S < 0, then (I4.14p would have k decay to zero no matter 
where it started; there would be no hint of a finite-amplitude instability in this case. 

Note that at the threshold of absolute stability we have V" = Q' c . In this case the 
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coefficient of k 5 in (14.141) is negative. In general, the calculation above then predicts that 
traffic will be in the regime where small- amplitude jams are seen, given by (14. 7p with k = 
(75/IOtK,' 2 ) 1 / 2 . This fails, however, if the threshold of absolute stability is at or near an 
inflection point of V s , as is the case in the OVM, the FVDM, and the COVM. 

V. NEAR AN INFLECTION POINT 

The picture changes dramatically if, in addition to being near the threshold of instability, 
the traffic spacing A is near an inflection point of the steady-state velocity function V s . 
Suppose the inflection point occurs at A = Aj, and V,'(Aj) is within order e 2 of the instability 
threshold, 

V:(A l ) = n c (A t ) + e 2 6 t . (5.1) 
Let A itself be within order e 2 of the inflection point, say 

A = A, + e 2 f3 

It follows that V^"(A) is also small, specifically 

K"(A) = eV/"^)/?. (5.3) 

Note that if V,' has a maximum at Aj, then V"'(A) is negative. Since K,"(Aj) vanishes, it 
follows that Vs'(A) differs from V^'(Aj) only by order e 4 . From (I3.7aj) we then find 

1 4- 2A 

D = ^— e 2 5 with 5 = 5 t - tt' c [3 (5.4) 

with corrections of order e 4 . As before, steady flow is linearly unstable for positive 5, linearly 
stable for negative 5. 

From (13. 6 p we see that if both D and V" are of order e 2 , then fx is of order e, and 
the leading-order terms are the three order-e terms on the right side that do not contain 
T-derivatives. It is appropriate, then, to return to the original ansatz (13.11) for x n , but retain 
the new definition (14.21) of the slow time T. As in the preceding section, we may substitute 
all this into the evolution equation (II .lft and expand (to fifth order in e), or we may simply 
take the general expansion (13. 6p and multiply each T-derivative by an additional factor e. 
After substituting (15 .4p for D and (15.31) for V" and canceling an overall factor of e, this latter 
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(5.2) 



procedure yields, to first order, 

f T = C u f zzz + Cyzfzfzz + C13/I + e 



2A 



f z J + ^V'"f z + CutAt + C 12T f z f T + C 21 f z 



2 2 

+ C 22a f z f zzz + C 2 2bfl z + C 2 zflfzz + C 2i f z 



(5 



To simplify the corrections, we again substitute the leading-order terms for each /y in the 
correction. We then differentiate with respect to z and write the resulting equation in terms 
of g to obtain 



9t 



1 + 3A 
6 



V',9* 



+ e 



1 + 2A P . yrl/ , 2 (1 + 2A) 1 



-C 22a gg zz + C 22b g z + C 23 g 2 g z + 7^K"V 



with the new coefficients given by 

C 22a = C 22a + C12TC11 + CixrC\2 = 

(1 + 2A) 2 



(1 + 2A) 2 /1 + A 



C 22 b 

C 2 3 



22b 



C\\tC_ 



12 



4 



4A 



C*23 + 3ChtCi3 + C12TC12 



1 + 2A 
4^ 



2-01 



(5.6) 

(5.7a) 
(5.7b) 
(5.7c) 



Note that both C 22a and C 2 2& vanish under the assumptions underlying the OVM and FVDM. 

The leading order of (I5.6P is the mKdV equation only if Q' c vanishes or is at most of 
order e. Thus, the mKdV equation only appears when both V" and Q' c are small for a given 
A. That is, the mKdV equation is not a generic feature of car-following models. Rather, it 
appears only in those cases for which the threshold of absolute stability occurs at or near an 
inflection point of the steady-state velocity function. If the inflection point and the threshold 
of absolute stability are not close together, then we obtain a KdV equation near the former 
(as we found in the preceding section) and the more general equation (15.61) near the latter. 
Equation (15.61) does, however, share some important features of the mKdV equation. 

In those cases where the mKdV equation does apply, its hyperbolic-tangent "kink" so- 
lutions are of particular interest. In the context of traffic, a one-kink solution describes a 
"domain wall" separating regions of uniform traffic flow with different speeds and densities. 
The mKdV equation has a one-parameter continuous family of such solutions, and one can 
carry out a solvability calculation to determine which are preserved by the correction terms. 
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Like the mKdV equation, the full equation (15. 6p also admits a one-parameter family of kink 
solutions. Remarkably, in the general case it is also possible to carry out the perturba- 
tion calculation to find members of the family that persist when the order-e corrections are 
included. 

One-kink solutions of the leading order of ( 15. 6 p are given by 

g (0 \z, T; k) = Mk tanh(kz + k 3 uT) , (5.8) 

where M and u are given by 

V'"M 2 - 3(1 + 2X)Q'M + 2(1 + 3X)V' = 0, u = -V'"M 2 . (5.9) 

6 

Since V" is negative, there are two solutions for M, one positive and one negative, and - u 
being negative - both move relative to traffic slightly more slowly than V^. As was the case 
for the KdV solitons, the various powers of k in (I5.8P reflect the invariance of the leading- 
order evolution equation under the simultaneous rescalings z — > kz, T — > k 3 T, g — > k^g, 
which in turn reflects the arbitrariness of the small parameter e. 

To see the effect of the order-e corrections on the kink solutions, we would like to carry 
out a perturbation calculation analogous to the one in the preceding section. That is, we 
write the first-order solution as 



9 



z,T) = g(°\z,T;k(eT)) + egM(z,T), (5.10) 



allowing the parameter k to evolve slowly in time, and substitute into (15. 6p to get 



~dk~ k 



1 + 2X a^5 + ^V"'q^ 2 - (1 + 2A) Vy°) 

2 i>z 1 2 3 y g sUzzz 



+ WW + C 22b {gf ] ? + C 23 g^g^ + 1^ (0)4 



(5.11) 



The linear operator C is now 



L ~dT 6 Vs dz* 2 Uc dz* 9 2 K dz 9 • [bAZ) 



Next we seek the function s(z,T) which is annihilated by the adjoint operator 

o = cu = -s T + ^v; Szzz - ^^W 0) ^ + \v:"g^ 2 s z . 



(5.13) 



The equation can, in fact, be solved using standard techniques; we find that s is given by 

s = s(kz + k 3 uT) with = sec h 2+p y, (5.14) 

dy 
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the exponent p being given by 



3M(1 + 2X)Q' C 



P ~ (1 + 3A)K' ' { } 

Note that for Q! c = we have simply p = 0, so s reduces to as in the preceding section. 

We now take the inner product of the function s with both sides of (15.111) to obtain the 
evolution equation for the parameter k, 



dg 



(0) 



dk 



k 



s, 



^gi% + ?K"V 0)2 - (1 + 2A)2 K' q [°) 



+ C 22a g^gf] + C 22b {gf) 2 + C 23 g^ + ^K"V° )4 



.(5.16) 



On the right side, we can change variables to y = kz + k 3 uT, integrate by parts (noticing 
that the boundary term does not vanish, because neither s nor the g^ 2 and g^ A terms in 
the square bracket vanish at infinity), and evaluate the resulting integrals explicitly in terms 
of gamma functions. Unfortunately the integral on the left side of the equation diverges, 
whether or not Q' c vanishes. Setting the right side of the equation to zero - the solvability 
condition - thus identifies the single k value for which the g^ given by (15.81) is the e — > 
limit of the perturbed solution. Specifically, this yields 

k 2 



[1 + 2X)5 + V" M(3 



5 + 2p 



;i + 2A) 2 (l+p)K' + 4 



(2+p)C 22b -C. 



22a 



M + 2C 23 M 



-(3 + p)K'"M 3 
o 



(5.17) 



For the FVDM, p, C 22a , and C 22b all vanish, and C 23 is equal to — (1 + 2A)V s "'/4. The result 
above then agrees with the result given by Ou et al. Jq|, provided we set f3 and V"" to zero, 
as they do (explicitly for and implicitly for V""). 

For k values which do not satisfy the solvability equation, the above perturbation calcu- 
lation is unable to determine what effect the correction terms have on the solution. Thus 
the solvability approach is silent on the question of how, and indeed whether, the "selected" 
kink is established. 



VI. DISCUSSION 

We have examined a very general car- following model embodied in Equation HI. ip . in 
which each car's acceleration is a general function of its current speed, the headway between 
it and the next car ahead of it, and the rate of change of the headway, subject to monotonicity 
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assumptions that rule out unreasonable descriptions of driver behavior and an assumption 
that sufficiently many derivatives exist. The most popular previously published models in 
this class, those for which the more delicate nonlinear calculations described below have been 
carried out, have assumed that the dependence of acceleration on the velocities is linear, 
with coefficients that are independent of headway 



Most car- following models are based on an "optimal velocity function" V S (A), written 
explicitly into the model, which gives the speed of uniform steady-state traffic flow in terms 
of the headway A between cars. In those models, uniform steady flow is found to be linearly 
unstable if the steady-state traffic speed is too sensitive to headway, i.e. if the derivative 
Vg'(A) exceeds some limit Q c which can be calculated in terms of the parameters of the 
model. It is also possible to have "absolute stability," which is the situation where V,'(A) 
remains below the instability limit for all A. The threshold of absolute stability is reached 
when parameter values are such that the maximum value of V^'(A) matches the instability 
limit Q c ; in the most popular models fl c is a constant, so this can only occur at an inflection 
point of the optimal velocity function. In the more general setting, the optimal velocity 
function must be defined implicitly from the acceleration function, as in (12.21) . The linear 
stability results for uniform, steady flow then continue to hold in the general case, with 
the crucial difference that the instability limit Q c can itself be a nontrivial function of the 
headway A. This implies, in particular, that the threshold of absolute stability does not 
generally coincide with an inflection point of the optimal velocity function. 

For steady traffic with a spacing near the onset of instability, linear stability analysis 
reveals that the perturbations that are unstable (or least stable) have long wavelengths. 
Using this observation, one often reduces a car-following model in this regime to a con- 
tinuum model, in which the spacing between cars is governed by a Burgers equation plus 
small corrections. In more restricted parameter ranges, different evolution equations for the 
traffic flow result: the Korteweg-deVries (KdV) equation (plus corrections) near the onset 
of instability, the modified Korteweg-deVries (mKdV) equation (plus corrections) near the 
threshold of absolute stability. The KdV equation admits a one-parameter family of soliton 
solutions, and the mKdV equation a one-parameter family of kink solutions; conventionally 
one carries out solvability analyses that determine unique soliton and kink solutions from 
among these families that are preserved in the presence of the correction terms. This pro- 
gram of reduction to nonlinear evolution equations has been carried out for a number of 
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tice models 
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car-followin g m odels |4|, la, |9|, |ll|, |12|, |15| that are special cases of (11.11) . as well as for lat 
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22J-I2JI, look-ahead 



2l| and models which incorporate finite time delays 
26] , look-back [18(, and a dependence on the acceleration of the next car [27]. All in 
all, it is not surprising that the derivations should continue to hold with those extra effects 
included, since the starting point of each derivation is the assumption that the traffic flow 
evolves on length scales long compared to the headway A and time scales long compared 
with the time scale on which an individual driver adjusts speed. 

We find that the derivation of the Burgers equation goes through at the level of generality 
of (II. ip . and that the coefficients in the resulting equation match those that have been 
obtained in the more specific models in the past. In particular, the coefficient of the nonlinear 
term in the Burgers equation is simply given by V^'(A). 

Near the onset of instability of steady flow, we find that the derivation of the Korteweg- 
deVries equation also goes through as in the more specific models, but if the instability limit 
Q c is a nontrivial function of headway A, then there is a crucial change to the coefficient 
of one of the correction terms which can lead to qualitative changes in the nature of the 
solutions. Specifically, if f^(A) vanishes, as in the previous models, then a solvability 
analysis of the one-soliton solutions of the KdV equation shows that there is a "selected" 
soliton only when uniform, steady flow is linearly stable. However, when that derivative is 
nonzero - and close enough to V^'(A) - then the opposite may be the case: a selected soliton 
exists only when the uniform flow is linearly unstable. 

Rather than using the usual solvability condition, we prefer to use a multiple-time-scales 



approach 
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14j to investigate the effect of the correction terms on the family of one-soliton 
solutions of the KdV equation. The advantage of this approach is that it not only identifies 
the "selected" soliton, but also shows what happens to the other, non-selected members of 
the family. Namely, the corrections drive a slow evolution of the soliton parameter, at a 
rate that vanishes for the selected soliton. The analysis shows that there are two possible 
behaviors. In one scenario, there is a selected soliton only when steady flow is linearly stable, 
and this selected soliton is unstable: smaller solitons decay to zero, larger ones grow. In this 
scenario, which is the one that obtains when fi' c (A) vanishes, the selected soliton marks the 
threshold of a finite- amplitude instability of the linearly stable steady traffic flow. Beyond 
this threshold amplitude, which goes to zero as parameters are adjusted toward the threshold 
of linear instability, jams grow large enough to invalidate the assumptions underlying the 
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derivation of the KdV equation. In the other scenario, a selected soliton exists only when 
steady traffic flow is linearly unstable, and the slow evolution of the soliton parameter that 
is caused by the correction terms drives the system toward this selected soliton. That is, 
there is a forward bifurcation from stable steady flow to a stable, small-amplitude jam. In 
either scenario, the evolution of the soliton is slow. Specifically, if the derivative of the 
optimal velocity function differs from its critical value by a small amount e, then the soliton 
adjustment takes place on a time scale of e~ 2 . Thus, even if the soliton is in the process of 
decaying to zero, it can appear to be a persistent feature of the traffic flow. 

When parameter values are such that steady flow is near the onset of instability and near 
an inflection point of the optimal velocity function, then the assumptions underlying the 
derivation of the KdV equation are not satisfied. In previous models this is what occurs 
near the threshold of absolute stability, and in this regime one reduces the model to a 
modified Korteweg-deVries equation with corrections (9), [lsj. We find, however, that this 
derivation does not go through in the general case. Rather, the result of the derivation is 
the mKdV equation plus an extra term. The existence of this term had been anticipated by 
Komatsu and Sasa [15|, and it had been seen by Hayakawa and Nakanishi [18| in a model 
including look-back; in our case the coefficient of that extra term is proportional to Q' C (A). 
Thus the mKdV equation only arises if this derivative is zero, or at least small, at the 
inflection point of V s - that is, if the threshold of absolute stability occurs near the inflection 
point. In the previous models, Q c was a constant, so this happened automatically; as we 
pointed out above, it is not the case generically. 

Like the mKdV equation, the new equation also admits a one-parameter family of kink 
solutions. Moreover, it is still possible to carry out a solvability analysis for the new equa- 
tion to determine which of these kink solutions are preserved when the correction terms 
are included. The multiple-time-scales approach, however, fails - as it does for the mKdV 
equation - because the slow time derivative of the kink parameter is multiplied by a co- 
efficient that turns out to diverge. This divergence, in turn, seems to come from the fact 
that changing the kink parameter involves changing the asymptotic values of the kink, i.e. 
changing the traffic spacing infinitely far from the location of the kink itself. This, and other 
aspects of the mKdV equation, will be the focus of future work. 

As noted above, the threshold of absolute stability occurs where the derivative of the 
optimal velocity function, V^'(A), matches the instability limit, Q c , at only one value of the 
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headway A, and is below Q c for all other headway values. We find that the behavior of 
traffic near the threshold of absolute stability can then be described by one of two reduced 
equations. If the threshold occurs near an inflection point of the optimal velocity function, 
that is if both V^"(A) and £l' c (A) are small, then the model can be reduced to a mKdV 
equation. This is the situation which has been explored in previous models. If, on the other 
hand, the threshold occurs away from an inflection point, then V S "(A) and f^(A) are close 
together, but neither is small. In this situation the KdV equation applies, and this is the 
regime in which the correction terms lead to a selected, small-amplitude soliton, representing 
a localized, small-amplitude "jam", when steady traffic flow is linearly unstable. 

In conclusion, we wish to point out that although the forms of the terms appearing in 
reduced equations like the KdV and mKdV are determined only by the scalings that go 
into their derivation, the coefficients of those terms are sensitive to details of the model. 
The behavior of the solutions of those reduced equations can be sensitive to the values, 
and especially to the signs, of those coefficients. Thus implicit assumptions made when 
formulating the model can constrain the coefficients (often constraining them to be zero!), 
and this can then have important quantitative and qualitative effects on the kind of behavior 
predicted by the resulting reduced equations. 
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